Thermalization of an impurity cloud in a Bose-Einstein condensate. 
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We study the thermalization dynamics of an impurity cloud inside a Bose-Einstein condensate at 
finite temperature, introducing a suitable Boltzmann equation. Some values of the temperature and 
of the initial impurity energy are considered. We find that, below the handau critical velocity, the 
macroscopic population of the initial impurity state reduces its depletion rate because of stimulation 
effects. For sufficiently high velocities the opposite effect occurs. For appropriate parameters the 
collisions cool the condensate. The maximum cooling per impurity atom is obtained with multiple 
collisions. 



In a Bose-Einstein condensate (BEC) many atoms pop- 
ulate the same quantum state. This system is, therefore, 
suitable to realize an atom laser That is obtained 
creating atoms in an untrapped or antitrapped level by 
means of a rf or Raman outcoupling. Elastic collisions 
between impurity atoms and the parent BEC have been 
observed at MIT . These collisions can deplete and 
heat the outcoupled atoms. In this work we study the 
thermalization dynamics of the impurity atoms, created 
by a rf or Raman outcoupling. We find the following 
main results: i) we evaluate the impurity thermalization 
and depletion time for some values of the temperature 
and of the initial impurity energy. The thermalization is 
not suppressed below the Landau critical velocity for a 
finite condensate temperature T, according to a previous 
work [j| ; ii) below the Landau critical velocity the macro- 
scopic population of the atom laser state slows down the 
impurity depletion, whereas, for sufficiently high veloci- 
ties the macroscopic population enhances the scattering 
rate, as observed in Ref. 0; iii) the collisions cool the 
condensate if the impurity initial energy e; n is below its 
thermal equilibrium energy, or even if this condition is 
not fulfilled, provided the thermal energy, divided by the 
chemical potential, is sufficiently small. These effects can 
be observed creating a succession of atom laser packets 
and observing a reduction of the impurity depletion rate 
in the last packets. This series can be used to reduce the 
condensate temperature, as proposed in Ref. In Ref. 
we evaluated the scattering rate of an impurity inside a 
condensate at finite temperature and found that for suit- 
able initial velocities the impurity acquired energy, on 
average, after the first collision. However processes with 
multiple collisions were not considered. Here we show 
that the maximum cooling per impurity atom is obtained, 
for T not too small, at thermal equilibrium, that is, with 
multiple collisions. It is important to understand that 
the first cooling is obtained by a weak rf or Raman out- 
coupling. In fact, the condensed fraction for a trapped 
condensate is n Q /n t = 1 - (T/T ) 3 [13, where n t , no are 
the number of overall and condensed atoms, respectively, 
and To is the transition temperature. The number of non- 
condensed atoms is nt — no = aT 3 , where a = n t /T^ is 
a constant . This means that the condensate temper- 
ature is given by the number of non-condensed atoms, 
with no regard to the number of the condensed ones. So, 



if the outcoupling does not create many impurities, the 
number of outcoupled non-condensed atoms is small and, 
consequently, the impurity temperature is lower than the 
parent condensate one. For this reason, it is possible to 
cool sympathetically a condensate creating a small im- 
purity cloud, as suggested in this work. 

To simplify the problem we consider an homogeneous 
condensate, in fact it is meaningful in a trapped system 
in the sense of the Thomas-Fermi approximation. If the 
number of outcoupled atoms is sufficiently low, the ini- 
tial thermal energy of the impurities can be neglected. 
We first introduce a Boltzmann equation for an impu- 
rity atom. A similar equation in momentum space is 
obtained, for a single particle and T = 0, in Ref. |6|. It 
is easy to show that the system goes toward the Boltz- 
mann equilibrium distribution. We extend our discussion 
to a cloud of impurity atoms introducing suitable bosonic 
stimulation terms. We obtain an equation of the same 
type of the quantum Boltzmann equation 0, B Hi] • It is 
verified that the Bose-Einstein distribution is the steady 
state solution. 

Let us consider an impurity atom in an uniform con- 
densate at temperature T. If the impurity energy is 
Ei, the probability that the particle scatters in a state 
with energy E f is T(E f ,Ei) = T 1 {E f ,E l ) + T 2 {E f ,E l ), 
where 
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hui® and TiLOq = 

flu® (hid® + 2fi) being, respectively, the initial impurity 

momentum, the condensate temperature in energy unit, 
the static structure factor and the energies of a free parti- 
cle and of a Bogoliubov quasiparticle of momentum q. M, 
fi and a are the atomic mass, the chemical potential and 
the scattering length for s-wave collisions between the im- 
purity atoms and the condensate ones, u is equal to nog, 
where no is the condensate density and g — 4:irh 2 ao/M, 
ao being the condensate s-wave scattering length. 

Ti and T2 are associated with dissipative and cool- 
ing processes, respectively. In the first case condensate 
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phonons are created, in the second one they are annihi- 
lated @ . With a little of algebra we find that T(E f ,Ei) = 
2('y/fi)G(2Ef/(j,,2Ei/fj), where 7 = Airn a 2 c, c = 
y/fi/M being the sound velocity, and 
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with (3 = (3[i/2 and Ae = ej — e^. W(£f, q) is 1 when > 
e f > 1 /e,/ , —1 when e/ > ej > 1/e/ and zero elsewhere. 
The Boltzmann equation is 



/" G(e,e)P(e)de- / G(e, e)P(e)de, 



(2) 



where P(2E/(i) is the density probability per unit energy 
and r = 7*. It is easy to demonstrate that P(e) oc ^/ee~^ € 
is a steady solution of Eq. J2J . 

The system dynamics depends exclusively on the prod- 
uct fifi — /i/(fcsT) and the initial probability distribu- 
tion. The condensate density and the scattering length 
are present in the rate factor 7 and, therefore, they in- 
fluence only the time scale of the dynamics. 

In Fig. ^1 we report the average energy acquired by 
the impurity and lost by the condensate thermal atoms, 
in k B T unit, i.e. (3[(e) - e in ] = 0[{E) - E m \. We have 
considered 1/(3 = 10 and some values (> 1) of the initial 
impurity energy. We can note that (E) = (3/2)fcsT for 
f — * 00, that is the average energy of a free particle at 
the equilibrium. In contrast, if we neglect the collective 
modes, a thermal particle of the condensate has an av- 
erage energy equal to ~ 0.514(3/2)fesT, that is, about 
one half of the classical value. It is easy to demonstrate 
that the number of lost thermal atoms SN, when the 
condensate energy decreases by 5E, is ~ 0.785E / (hsT). 
Therefore, if the impurity acquires an energy of ~ 1.3k bT 
a thermal atom jumps into the condensate state, bal- 
ancing the loss due to the impurity creation. A nearly 
complete thermalization occurs for t ~ 0.5, that is, for 
t t h — l/(2<7o™oc), where erg = 47ra 2 . For a 23 Na con- 
densate and a typical density of 10 14 ~ 10 15 cm -3 (T ~ 
350 ~ 3500nK for 1/0 = 10) we have t th ~ 10 ~ 0.3ms. 
We can see that the thermalization time is independent 
on the initial velocity. In fact, the impurity loses quickly 
the memory of its initial energy. In the inset of Fig.^ wc 
report the relative density Pi{r) = P(e irl , T)/P(e irl , 0) as 
a function of r. The initial state population decreases by 
60 -i- 70 per cent for r ~ 0.1, i.e. for t = 0.06 -7- 2ms. We 
note that the depletion rate is smaller for a higher initial 
velocity, contrary to the zero temperature case. In fact, 
for T = 0we expect that the scattering rate grows upon 
increasing the initial velocity. At finite temperatures, 
to reduce the laser mode depletion, it is convenient to 
have a large initial velocity, because the scattering rate 
is smaller and, furthermore, the impurity leaves the con- 
densate earlier. For a sufficiently high velocity we expect 
that the scattered impurity fraction is not enhanced by 
the finite temperature. 




FIG. 1: Average energy acquired by the impurity in fcsT 
unit as a function of r, for some values of e,: n . In the inset 
is reported the relative density Pi(r) as a function of r. (a) 
1/(3 = 10, i.e. k B T = 5/i. (b) 1/(3 = 2, i.e. k B T = fi. 




FIG. 2: Average energy acquired by the impurity in fj,/2 unit 
as a function of r, for some values of ei n and for 1/(3 — 0.6-r-l. 



If the condensate temperature is reduced, the scatter- 
ing rate decreases and approaches the zero temperature 
value. For T = only the spontaneous scattering con- 
tributes and when Ci„ < 1 the impurity thermalization 
is absent because of the condensate superfluidity. For 
ksT = [i (/3 = 2) the temperature can be still important. 
In fact, we can see in Fig.^ 3 that a nearly complete ther- 
malization occurs for r ~ 2, that, for hq — 10 15 cm~ 3 , 
corresponds to t ~ 1.2ms. The initial velocity is be- 
low the Landau value, so the first impurity scattering is 
due to the finite temperature. The function G(e/,ej) is 
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equal to zero for ej < 1/ej, therefore, if the initial en- 
ergy is too low, the particle jumps in a state with a large 
energy. However the probability rate of this process is 
small, because of the presence of the exponential in G. 
So we have a slowing down of the thermalization rate 
(see Fig. \iy>, solid line) when the initial velocity is very 
low. For ein > 0.4 the cut-off of the constraint Ef > 1/ti 
has a minor effect and the thermalization time becomes 
again independent on the initial velocity. Also the de- 
pletion rate of the relative density Pj(r) is reduced by 
the cut-off for e in = 0.2. Increasing e in the depletion 
rate grows, however for ei n > 0.6 it slows down again, 
as it occurs for j3 = 10 (see Fig. ^). For n in = 0.4 



and no = 10 
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15, 



,-3 



the initial impurity velocity 



is ~ 0.3 4- lcm/ s, therefore, the average distance before 
the first collision (r ~ 1) is about 4 • (10 -2 4 10~ 3 )mm. 
The radius of a condensate with the same central density 
is ~ 0.5 4 2 • 10 _1 TOm when the trap frequency is 20Hz. 
To contrast the gravitational acceleration of the impu- 
rities, both the parent and outcoupled condensates can 
be trapped in the radial direction with a far off-resonant 
Gaussian laser beam. In the axial horizontal direction a 
magnetic field can be used to trap only the parent con- 
densate. We can set the trap in a such way to permit 
that the impurities leave the condensate after the ther- 
malization time. 

In Fig. |2 we report the same curves of Fig.^ but the 
energy is in /i/2 unit and 1//3 = 0.6 4 1. In this case 
the energy reaches a maximum at an intermediate time. 
For 1/(3 = 1 and = 0.6 4 1.0 the maximum is reached 
at t = 10 4 20, that, for no — 10 15 cm -3 , corresponds 
to t = 6 4 10ms. We find that the evolution of the 
initial state population is nearly independent on ej„. At 
t = 10 the initial state is depleted by 40%. To cool the 
condensate it is not necessary to reach the maximum, 
but for smaller times we obtain a minor reduction of the 
condensate temperature per atom. Note that for 1/(3 = 
0.6 the equilibrium thermal energy is 0.9, in p/2 unit. 
Therefore, for ej„ = 0.9 4 0.95 the initial energy is equal 
or higher than the equilibrium value. Even so, at an 
intermediate time the impurity acquires energy and the 
condensate cools. That occurs because when the velocity 
is below the critical value the first scattering can only 
enhance the impurity energy. 

Let us suppose now that there are many impurities, 
then many atoms can be in the same energy level and, 
therefore, it is necessary to consider the stimulation ef- 
fects. If we suppose that the initial impurity velocity 
is isotropic, we can make the following modification in 
Eq. © 



dpje) 
air 



\p(e)/p(e) + l] / G(e,e)p(e)de 



[p(e)Me) + l]G(e, e)p(e)de = B(p), 



(3) 



where p(e) = P(e)/V, V being the system volume, and 
p(2E/p) = {2m) 3 /' 2 h~ 3 (2Tr)- 2 VE is the level density per 
unit volume jToL ITU ] . It is easy to demonstrate that the 




FIG. 3: Relative density Pi(r) as a function of r for 1//3 = 2 
and some values of ni/n t h, with (a) ei n = 0.2 4 0.4, (b) = 
1.5 and e™ = 1 in the inset. 



Bose-Einstein distribution p(e) oc ^/e{e^^^ — 1) _1 is a 
steady solution of Eq. © , p, being the normalized chem- 
ical potential of the impurity cloud. 

We isolate in p(e) the contribution of the initial popu- 
lated state with energy e in . The associated distribution, 
that is very narrow, has a width equal to Se = 2/{V p\x). 
We indicate with its height. 

Eq. © becomes 



Me) 
dr 



dpi 
dr 



B(p)+8e[p{e)/p{e) + l]G(e,e ln )p l 

-Se[pi/p(e in ) + l]G(e in , e)p(e) (4) 
[Pi/p(ei) + !] / G(e in ,e)p(e)de 



\p(e)/p(e) + l]G(e,e in ) Pi de. (5) 



The quantities 5eG(ei n , e)p(e) and / G(ei n , e)p(e)de in the 
first and second equation, respectively, can be neglected. 
Therefore Eqs. (|4I5() become 



Me) 
dr 



B(p) + -{\p{e)/p(e) + l]G{e,e in ) 
P 

-G(e in ,e)p(e)/p(e in )}ni (6) 

= -7»(p)»i) ( 7 ) 
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FIG. 4: Same as with Fig.Q but for an impurity cloud. Some 
values of m are considered, for = 2 and e in = 0.6. 

where rij is the spatial density of the impurity cloud and 

liip) = J W)/P® + l]G(e,e in )de- 

l/p(e,n) J g(e ln ,e)p(e)de = J{G(e, e in ) - 

W(e in ,e)p(ei n )- 1 [l - (1 + Ae 2 /Ay^p^jde, (8) 

Ae being ej n — e. The first term of 7.; is the depletion rate 
of the initial state when the stimulated effects of the im- 
purities are negligible. The second term, that does not 
depend directly on the temperature, is proportional to 
the density p(e) and it is important when there are many 
impurities. This term is negative when e„ < 1, therefore 
in this case it reduces the depletion rate ji. In fact, the 
scattered atoms can have a high probability to scatter 
back into the original macroscopically populated state. 
On the contrary, when £j„ > 1 the second integrand is 
negative for e > e, n and positive for e < e;„; if e,-„ is suffi- 
ciently high the impurity stimulation effect enhances the 
depletion rate, as it has been experimentally observed 0- 
In our calculations we have considered the condensate as 
a thermal reservoir, however if we have many impurities 
the collisions can cool or warm the condensate, conse- 
quently the depletion rate can decrease or increase. Our 
approximation is reasonable when the density of outcou- 
pled atoms is small with respect to the density of thermal 
atoms n t h or, otherwise, for sufficiently small times, when 
the collided fraction is small. We have approximatively, 
if we neglect the condensate collective modes, 0] 

nth = (4 7 r 2 ^)- 1 m 3 / 2 ( M //3) 3 / 2 (9) 
In Fig.|3t we report the relative density Pj(r) for l/j3 = 



2, ei n — 0.2-^0.4 and some values of nil nth = ni(0) /nth- 
It is easy to demonstrate that these curves are indepen- 
dent on [i and, therefore, on a and no. In fact, with 
p(e) oc y/Ji and nj oc /i 3 / 2 , the Boltzmann equations are 
invariant with respect to modifications of the chemical 
potential. For no = 10 14 cm~ 3 and uq = 10 15 cm~ 3 we 
have nth — 1-1 ■ 10 12 cm~ 3 and nth — 3.5 • 10 13 cm~ 3 , re- 
spectively. We can see that for nj/n t h = 10 the depletion 
rate is considerably reduced by the stimulation effects. In 
Fig. 03 we plot the same curves but for e in = 1 (inset) 
and 1.5, respectively. These curves are reliable only at 
the initial times because of the cooling of the condensate. 
Furthermore, if collisions create many thermal impurity 
atoms, a lot of impurities could condensate into the state 
with zero energy, but our equations are not suitable to 
describe this phenomenon, because this state has to be 
treated separately from the energy integral of Eq. Jfjjl. 
If the number of scattered atoms is lower than n t h we 
can suppose that this condensation does not occur. Note 
that the depletion rate is lowered also for ti n — 1.5, that 
is above the Landau value. Instead, for Ci n = 2 we find 
that the depletion rate is enhanced. Increasing the tem- 
perature the threshold of e, n increases. 

Finally, we report in Fig. 0] the same curves of Fig.^ 
but accounting for the stimulation effects of the impurity 
cloud. We have considered some values of n/, for 1//3 = 2 
and — 0.6. The maximum acquired energy is reduced 
increasing n/. 

In conclusion, we have evaluated the thermalization 
time of impurity atoms inside a condensate at finite tem- 
perature and we have found that their thermalization, 
that occurs also below the Landau critical velocity, can 
cool the condensate. In fact, a small Raman or rf outcou- 
pling creates a condensate in another level with a lower 
temperature, that cools sympathetically the parent con- 
densate. We have shown that this effect can be studied 
experimentally. It could be used to remove the resid- 
ual thermal atoms. Furthermore, we have accounted for 
the stimulation effects of the impurity cloud and we find 
that below the Landau critical velocity the macroscopic 
population of the initial state can reduce drastically its 
depletion rate. For sufficiently high velocities the oppo- 
site effect occurs, as observed in Ref. 0- The threshold 
velocity depends on the condensate temperature. Our 
analysis is also important for the development of intense 
atom lasers. 
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